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ABSTRACT 



The purpose of this effort was to gain an insight into 
the Bi-phase cycle concept and to develop an experimental 
facility to test the key component of this power cycle — the 
Bi-phase nozzle. A system was developed and checked out that 
allows testing different Bi-phase nozzle shapes. The experi- 
mental system is provided with sufficient instrumentation to 
allow for the determination of the nozzle performance. In 
addition to the standard gas and liquid flow rates, measure- 
ment provisions exist for survey of the nozzle exit plane. 
Measurements can be made for the gas/liquid ratio and 
velocities at the nozzle exit. 
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I. 



INTRODUCTION 



Low speed prime mover systems in marine applications 
are attractive in that they negate requirements for heavy, 
voluminous and noisy reduction gearing. Direct drive con- 
ventional steam turbines and conventional electric drives 
have been tested and utilized with only limited success. 

What is required is a machine with low specific speed that 
is compatible at the same time with existing heat source 
systems. One option that may meet these requirements is 
the Bi-phase prime mover. The Bi-phase turbine is a 
machine having low specific speed and high torque. As 
such, it is ideally suited for marine propulsion. The Bi- 
phase turbine can drive the propellers directly without 
gear reduction. This feature leads to improvements in 
weight, volume and noise reduction. 

In addition, certain aspects of the Bi-phase concept 
may allow increased power production efficiencies. It may 
be possible to achieve near constant temperature expansion 
in the two-phase nozzle thus realizing a degree of 
"Carnotization" of the basic Rankine cycle. A detailed 
description of the Bi-phase power plant is given in Section III. 

A key element in the development of the Bi-phase prime 
mover is the Bi-phase nozzle. A Bi-phase nozzle is a flow 
accelerator in which fluid particles or droplets are 



8 



accelerated by a vapor. The liquid and vapor may be 
chemically identical or dissimilar. The expanding vapor 
imparts kinetic energy to the liquid droplets. During 
this expansion, the droplets may also transfer heat to or 
from the vapor, resulting in a near isothermal expansion. 

In the nozzle, most of the available thermal energy of the 
vapor is converted to kinetic energy of the liquid droplets. 
The nozzle is converging-diverging, a requirement resulting 
from the low sonic velocity in two-phase mixtures. Computer 
codes [Ref. 1] exist that yield nozzle exit velocities to 
reasonable accuracies provided fluid properties are known, 
liquid particle size can be controlled, and one dimension- 
ally in the flow field is assumed. 

Tangren, Dodge, and Siefert derived the equations for 
isentropic flow of an immiscible two-component mixture 
with equal temperatures of the liquid and the gas. The 
first extensive experiments on two-phase nozzles were 
performed by Elliott [Ref. 1], 

The purpose of this effort is to develop an experimental 
facility that will allow evaluation of two-phase nozzle 
concepts and associated liquid-vapor interaction phenomenon. 
The ultimate goal is to provide experimental capability to 
validate analytic models of Bi-phase nozzles. To this end 
an experimental system was designed that includes a high 
degree of flexibility to allow for rapid flow passage 
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geometry variation. The instrumentation installed has the 
capability to sample liquid content at different points in 
the nozzle exit plane. 
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I I . BACKGROUND 



Turbine type prime movers, whether based on the Rankine 
or Brayton cycles, have many advantages. Among these are: 
high torque at low speed, low vibration levels as compared 
to internal combustion energy converters, and in some cases 
very high power densities. Practical considerations result 
in prime mover power outputs which are at relatively high 
rotational speeds. 

Power demand devices such as propellers and pump jets 
have rotational speed requirements set by considerations 
considerably different than those of the prime mover. Un- 
fortunately the rotational speeds of prime mover and power 
user seldom match. As a consequence, a power conditioning 
device, a gear box or reduction gearing is invariably 
mandated. In certain installations the reduction gearing 
may weigh more than the prime mover. 

From a purely thermodynamic consideration one would 
want the maximum cycle temperature to be as high as possi- 
ble. However, where the working substance must be 
contained by physical structures such as combustors, 
boilers, nozzles, and turbine blades, the cycle temperature 
limit is set by material considerations. Furthermore, 
secondary aspects of a prime mover, i.e., safety, may 
further restrict the range of cycle parameters. A case in 
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point is the Navy nuclear power plant that does not admit 
superheat. The turbine inlet conditions are those of 
saturated vapor. This constrains the maximum cycle effi- 
ciency to comparatively low values. 

A relatively new prime mover concept that may have the 
capability to circumvent some of the above identified diffi- 
culties is the Bi-phase prime mover. The technology base 
for the Bi-phase concept is founded in earlier work at the 
Jet Propulsion Laboratory on liquid metal power generating 
concepts. 

The baseline Bi-phase concept is still the Rankine 
cycle (Fig. 1) . The differences are subtle but of extreme 
importance. The working fluid is vaporized in the heat 
addition component of the cycle. This may be a conventional 
boiler generating high temperature and pressure steam. 

Normally the working fluid would then enter a nozzle/turbine 
system and energy would be converted to a mechanical form. 

In the Bi-phase system, however, a second fluid is intro- 
duced at boiler temperature and pressure just prior to the 
point where the steam enters the nozzles. This second 
fluid may have a vaporization temperature higher than water 
and thus will remain in the liquid state throughout the 
cycle. In the nozzle proper the vapor constituent acceler- 
ates and in so doing drags the droplets of the liquid phase 
to higher velocity (Fig. 2) . This energy exchange of enthalpy 
drop of the vapor to kinetic energy of the liquid is the 
major departure from the conventional steam power plant. 
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BI-PHASE CYCLE 




Fig. 1. Temperature-Entropy Diagram 
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Fig. 2. Liquid Entrainment Photograph 
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Downstream of the nozzle the high kinetic energy liquid 
and the expanded vapor are separated. The vapor is condensed 
and returned to the boiler. The liquid is "parked" on a 
rotating wheel and energy is extracted by an impulse or pitot 
type turbine. After interaction with the hydraulic turbine, 
the fluid is returned for injection at the Bi-phase nozzle 
mixing section. 

It is the Bi-phase aspect of the system that is at the 
base of some of the potential advantages of the concept, 
which are: 

1. The energy imparted to the liquid phase is in a 
kinetic form of a high density mass. Hence, the 
liquid velocities involved are relatively low as 
compared to velocities of the vapor. Thus the 
output of the impulse hydraulic turbine will be at 
high torque and low rotational speeds — more in 
keeping with rotational speeds of driver elements. 

2 . The energy extraction from the steam and the energy 
addition to the liquid takes place in an inherently 
small device — the Bi-phase nozzle. Thus, the high 
temperature, pressure, and velocities are confined 
to a high energy density component. 

3. Perhaps the most important aspect is the phenomenon 
that takes place in the nozzle proper between the 
liquid and vapor and how this interaction impacts 
the overall cycle performance. While it is normal 
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practice to assume "isentropic" expansion in a 
nozzle, it is not likely that this is what happens 
in a Bi-phase expansion. In all probability thermal 
energy is exchanged between the liquid and vapor 
constituents. This, in fact, may be the most 
attractive aspect of the Bi -phase concept. If 
thermal energy can be added to the vapor phase 
from the liquid, then the expansion can be made to 
approach an "isothermal” process. The net result 
is a degree of "Carnotization" of the basic Rankine 
cycle that can only lead to increased overall cycle 
efficiencies . 

Application of two-phase turbines to Navy propulsion has 
been examined in two design studies. References 2 and 3. 

In both cases, the following potential advantages were found 
when the Bi-phase turbine was compared to advanced steam 
turbine systems. 

1. High efficiency — Full load output power performance 
gains ranging from 20 to 50 percent were found. 

2. Direct drive — Direct drive at speeds ranging from 
90-4500 rpm was found possible with a single-stage 
turbine . 

3. Reduced volume — Volume reductions of 30 percent 
were estimated. 

4. High part-load efficiency — Variable mass ratio 
enabled part-load (cruise) efficiency gains of as 
much as 100 percent. 



16 



An independent Navy study of a large Bi-phase marine 
turbine arrived at these same conclusions. These potential 
advantages lead to many possible applications for marine 
use. Most prominent among these advantages are those 
applications where space and weight are the most important 
considerations, for example, torpedoes, submarines, SES, 
SWATH, etc. 

It is thus apparent that the key element of a Bi-phase 
power system is the Bi-phase nozzle. It is this component 
that is the subject of this thesis. Specifically, the goal 
is to develop a research capability to study the flowfield 
and energy interaction in a Bi-phase nozzle. 
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III. BI-PHASE CYCLE THERMODYNAMICS 



The Bi-phase cycle is basically a variation of the con- 
ventional Rankine steam plant. In reality, the "Bi-phase 
aspect" should be a multicomponent cycle as there is no 
requirement that the components (vapor and liquid) of the 
working substance be chemically identical. 

In a conventional Rankine cycle (Fig. 1) the liquid 
phase is introduced to the heat additive component of the 
system at the maximum cycle pressure. Energy is added and 
the substance is vaporized. This is a constant pressure 
energy addition process and follows a constant pressure line 
on the T-S coordinate system. The maximum cycle temperature 
may be well into the superheat region. Expansion, often 
idealized as isentropic, follows. During this pressure drop 
from boiler to condenser pressures, work is extracted by a 
turbine via a volume expansion. Constant pressure heat 
rejection returns the working substance to a liquid phase 
which then is pumped up to boiler pressure and reintroduced 
to the heat addition component. 

A common way of depicting the maximum cycle efficiency 
for any cycle is to apply the Carnot efficiency criteria. 
This Carnot efficiency depends only on temperatures of heat 
addition and heat rejection. 
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Where n n = Carnot efficiency 

= Condenser temperature (°K) 

T = Boiler temperature (°K) 

Thus, a cycle operating between source-sink temperatures of 
232°C and 64°C respectively has a maximum energy conversion 
of 33%. The efficiency for a Rankine cycle with the same 
temperature limits and saturated vapor turbine inlet condi- 
tions yields an ideal efficiency of 29%. This efficiency 
is based strictly on thermodynamic considerations and does 
not consider irreversibilities such as friction losses, 
temperature drops in heat exchangers, and hydrodynamic losses 
in the cycle mechanical components. Furthermore, from 
practical considerations, the last stages of the expander 
see a very wet working fluid = 75% quality. This results 
in erosion and impact losses leading to further decrease 
in the power generator effectiveness. 

It is the Bi-phase cycle variation that has the potential 
to circumvent the above identified difficulties of the 
Rankine plant. 

The Bi-phase cycle adhieves improved efficiencies because 
the process involved allows the development of a near Carnot 
cycle. Thus, the Bi-phase concept achieves a degree of 
Carnotization of the Rankine cycle. This Carnotization 
takes place in the Bi-phase nozzle process. While in a 
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Rankine plant the expansion process is idealized as isen- 
tropic; in the Bi-phase system, it can be caused to approach 
an isothermal process. 

There are basically two Bi-phase techniques whereby a 
Rankine cycle can be caused to approach the Carnot cycle. 

A. TWO COMPONENT SYSTEM 

One method to Carnotize a steam cycle is to add a second, 
chemically different component that remains a liquid through- 
out the cycle process. This component provides a source of 
heat for the expanding vapor producing a near-isothermal 
expansion. In Fig. 1 it can be seen that the temperature at 
the end of the. nozzle expansion is 218°C compared to 64°C 
which would occur in the Rankine cycle. This results in 
a cycle efficiency of 31% versus 29% for the Rankine cycle. 
This cycle was invented by D. G. Elliott of Cal Tech for 
liquid metal MHD . [Ref. 1] 

The two-component cycle uses a low vapor pressure liquid 
and a high vapor pressure liquid having different chemical 
makeup. Some fluid combinations which have been considered 
are steam-Krytox, steam-Caloria, steam-lead bismuth eutectic, 
and Dow-therm-Therminol . The basic advantages of a two- 
component engine are low rpm, high efficiency with no steam 
extraction and a more compact turbine. Disadvantages 
include larger heat exchangers and the complexity of handling 
two fluids. Figs. 3 and 4 depict two possible methods of 
utilizing the two-component approach. 
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Fig. 3. U-Tube Turbine Engine 
21 





Fig. 4. Impulse Turbine Engine 
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D. SINGLE COMPONENT SYSTEM 

In this cycle the expansion process is ideally paralled 
to the saturated liquid line that outlines the left hand 
side of the enthalpy dome. In order to achieve this type 
of Carnotizing, heat has to be removed during expansion. 
While it may be possible to do that through the walls of 
the nozzle, a more conventional way is to extract steam in 
steps. That means that a scheme of pressure staging is 
required, where after each expansion stage steam is ex- 
tracted and used to preheat feed water by condensing all of 
the extracted steam. 

The single component system has the advantage of elimin 
ating oil or heat transfer fluid as the second component. 
The two-phase mixture in this case consists of steam and 
water. Typical volume ratios of water to steam are 10:1 
or less so this cycle produces somewhat higher values of 
turbine speed. 

In order to obtain high efficiency several expansion 
stages are used. The five stage system resulting in the 
temperature-entropy variation of Fig. 5 is shown schematic- 
ally in Fig. 6. Steam from the heat exchanger flows to the 
first stage nozzle. A portion of the steam is also used 
to heat the water flowing to the first stage nozzle. The 
mixture of primary steam and water is expanded in the 
nozzle from an inlet vapor quality of 15% to an exit 
quality of 13%. The two-phase mixture drives a rotary 



23 



400 r 



300 



TEMP 

(°C) 

200 . 



100 



j i l I I 1 I L 




ENTROPY (kj/kg) 

Fig. 5. Single Component Bi-phase Cycle 



24 



VO 



p^l 



VO 




VO 



u 

<D 

cn 

c 

(U 



h 


G 

0 






s'Lr — — ■ ■ 


u 

in 


x ^ — — 





lO 




ro 



CN 



e 

(d 

(U 

-p 

in 



Fig. 6. Wet Steam Five Stage Regenerative Turbine Engine 
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Pump 



separator turbine. Part of the separated steam is extracted 
and flows to a feed water heater. The remainder of the 
steam is remixed with the water exiting the first stage 
turbine. This mixture is expanded in the second stage 
nozzle from an inlet vapor quality of 12% to an exit 
quality of 18%. 

The two-phase mixture drives the second stage rotary 
separator turbine, where additional steam extraction occurs. 
In the final stage, steam exhausting from the turbine flows 
to the condenser. The condensate is remixed with the 
liquid exiting the last stage turbine. The feed water 
flows through several stages of a feed water heater where 
the extraction steam is condensed on the feed water. 

Five stages are shown in Fig. 6. The actual number of 
stages will be determined by a tradeoff between efficiency 
and cost and complexity. 
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IV. BI-PHASE NOZZLE THEORY 



The Bi-phase nozzle is the key element in the Bi-phase 
engine concept. While initially it may appear that the 
design of the nozzle profile is quite straight forward, 
closer examination will reveal that the question is con- 
siderably more complex. The complexity steams from the 
nature of the working fluid itself. It is a mixture of 
liquid particles and a compressible gas or vapor. 

The complicating factors can be grouped into three 
categories : 

1. The working fluid has two constituents significantly 
different in density. Thus, there is an aerodynamic 
interaction between the components (liquid droplets 
are accelerated by the vapor) . 

2. The shape and size of the liquid particles is 
governed by the forces on them. In this particular 
case the surface tension and inertia forces are of 
key interest. 

3. Thermal interaction between the working fluid 
components can range from none at all to phase 
change of the liquid to its vapor form. 

Another factor of interest in Bi-phase flow is the 
question of the velocity of propagation of a disturbance, 
i.e., velocity of sound. While the velocity in a particular 
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vapor and in a particular liquid are well established, the 
situation in a mixture of vapor and liquid is considerably 
different . 

Consider the following discussion from Olson [Ref. 4], 
The acoustic velocity in any medium is c = /K/p , where K is 
the modulus of elasticity for the medium and o is the 
density. For gas-liquid mixtures the acoustic velocity 
becomes less than that for either the liquid or gas alone. 
For a liquid with a small concentration of gas nuclei, the 
elastic modulus of the mixture is reduced, with no appre- 
ciable reduction in density, and thus the acoustic velocity 
is reduced. For a gas with minute liquid droplets, the 
density of the mixture is increased with no appreciable 
change in elastic modulus, and again the acoustic velocity 
for the mixture is reduced. The velocity of sound in a 
14 percent water and 86 percent air mixture is about 30 m/s 
(depending on the pressure, temperature, and impressed 
frequency) as compared with about 1460 m/s for water and 
about 335 m/s for air. 

The flow phenomenon of a Bi-phase mixture has been 
analyzed in Ref. 1. It is repeated as follows. The 
problem is illustrated in Fig. 7. A spatially uniform two- 
component mixture of liquid drops and gas enters a nozzle 
at high pressure and low velocity and expands to low 
pressure and high velocity. The objective of the analysis 
is to determine, for a specified pressure, the drop diameter 
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D and the temperatures T and T_ , velocities V and V T and 

g L g L 

flow rates m^ and m^ of the gas and liquid phases, respec- 
tively, at each station in the nozzle given the initial 

values of D, T , T t , V , V T , the total flow rate, and the 
g L g L 

properties of the fluids. 

The five relations employed to compute the five unknowns 

D, Tg, T l , V , and V^, are (1) the momentum equation for the 

mixture, (2) the energy equation for the mixture, (3) the 

drop drag equation, (4) the drop heat transfer equation, and 

(5) the drop breakup criterion. Solubility and vapor 

pressure relations provide the flow rate ratio m /m . 

g l 

A. ASSUMPTIONS 

The assumptions employed in the two-component analysis 
are as follows: 

1. The liquid is uniformly dispersed as spherical drops 
all of the same diameter. 

2. The drops break up to limit the Weber number to 6. 

3. There are no external forces acting on the two-phase 
mixture other than pressure and wall shear, and 
there is no heat transfer to or from the mixture. 

4. The flow is one-dimensional. 

5 . The drops are large enough for the surface curvature 
to have negligible effect on the vapor pressure of 

the liquid and for the surface energy to be negligible. 

6. The drops are isothermal. 

7. The gas mixture obeys the additive-pressure law. 
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The partial pressure of the predominantly liquid 
component in the liquid is given by Raoult's Law. 

9. The concentration of the predominantly gaseous 
component in the liquid is given by Henry’s Law. 

10. The volume of the liquid solution is equal to the 
sum of the volumes of the pure liquids. 

Assumption 1 restricts the analysis to nozzles having 
spatially uniform injection of the liquid into the gas and 
operating at gas-to-liquid volume ratios greater than 
unity. Assumption 2, the drop breakup criterion, states 
that drop diameter is limited to a value D for which W 



W = p V D/2a = 6. Thus 
e w g g 



max 



12a 
P V 2 

g g 



( 1 ) 



where is the gas density, is the slip velocity V -V^, 
and a is the liquid surface tension. The form of Eq. (1) is 
physically reasonable in that the Weber number is proportional 
to the ratio of stagnation pressure PgVg/2 to surface 
tension pressure 4a/D. Hence, a drop would be expected to 
flatten and break up at a sufficiently high value of We. 

This has been verified experimentally and the critical Weber 
number found to be 6, within a factor of about two. An 
additional restriction is that for actual breakup to occur, 
the time spent at a Weber number exceeding 6 must be longer 
than the natural period of oscillation of the drop, 
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3 ^ 

tt(p t D /a) 2 /4, where p is the density of the liquid. This 
requirement is met only in two-phase nozzles longer than 
about 10 in. and Assumption 2 may cause the analysis to 
overestimate the exit velocity by increasing amounts as 
the nozzle length decreases below 10 in. 

Assumption 3 excludes magnetohydrodynamic and mechanical 
body forces. The exclusion of wall heat transfer is correct 
for the insulated nozzles of interest for power systems. 

In addition, the relatively high velocity results in short 
residence lines in the nozzle proper. 

Assumption 4 is closely met in practical nozzles since 
good performance requires small wall angles, large throat 
radius of curvature, and uniformly distributed injection 
of the fluids at the nozzle entrance. 

Assumption 5 is valid for the drop sizes of 0.001 to 
0.010 in. produced by the Eq. (1) breakup criterion. 
Assumption 6 is valid because of the rapid internal circula- 
tion in drops. Assumption 7 introduces negligible error 
in most cases of practical interest since the vapor pressure 
of the liquid is small and needs only to be evaluated 
approximately . 

Assumptions 8, 9, and 10 are either valid, or cause 
little error, for fluids of low miscibility, which are the 
fluids of interest. 
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B. DERIVATION OF EQUATIONS FOR FREE-STREAM FLOW 



1 . Continuity 

Referring to Fig. 1, the nozzle flow area A is equal 

to the gas flow area m /p V plus the liquid flow area 

y y y 

m L /p L V L . Thus, 





( 2 ) 



where r is the mass mixture ratio m^/m^. 

2 . Momentum 

By Assumption 3, the only force acting on the free- 
stream flow is that due to the pressure gradient. If M is 
the momentum flux at flow area A, the change in momentum 
flux across pressure increment dp is 



dM = 



-Adp 



(3) 



The momentum flux can be written as the sum of the 
momentum fluxes of the gas and liquid. Thus, 



M 



m V 

g g 



+ m lVl 



(4) 



If the flow were allowed to continue at constant 
pressure, and would become equal to each other at the 
mass-weighted mean velocity V. Since for this process, 
dM = 0, the value of V is given by 




m, ) V = mV 

i g g 



m i v i 



(5) 
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or 



V + rV, 

V = ? - 1 

1 + r 



( 6 ) 



Thus, the momentum flux can be written 



M = (m + m. ) V 

g 1 



(7) 



Since m g + m^ is constant, the change in momentum flux is 



dM = (m + m. ) dV 

g 1 



( 8 ) 



Substituting Eqs. (3) and (2) into Eg. (3) , dV can be written 



dV = 



1 + r 



1 + r P V p.V 



. g g "ii. 

The slip ratio is defined as 



dp 



(9) 



s = Vg /V = (V - V ] _)/V 



( 10 ) 



This equation can be combined with Eq. 



(6) to give V and 



V L in terms of V; 



v * n 1 + 

V 1 ‘ i 1 - T*r) 



V = aV 



V = bV 



( 11 ) 

( 12 ) 



The gas density can be expressed as 



p = W p/RT 

g g g 



(13) 



where W g is the effective molecular weight of the gas mixture 
and R is the universal gas constant. Eq. (13) is the 
definition of the effective molecular weight W g , which is the 
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quantity that gives the actual gas density when substituted 
in Eq. (13) . 



Substituting Eqs. (11)-(13) into Eq. (9), the 
differential momentum equation is 



2VdV = dV 2 



( RT 

[ SL + 

1 + r l aW 

\ 9 




(14) 



The quantities a and b are slowly varying because s is 
typically only 0.1 to 0.3 and slowly varying. The quanti- 
ties r, Tg, W and p L are also slowly varying. Integrating 
Eq. (14) over a pressure increment Ap, for which a, b, r, 

T , W , and p L are constant to within the desired accuracy, 

. 2 . 

the change mV is 




(15) 



All quantities other than pressure can be taken outside the 
integral and evaluated at their mean values (denoted by 
subscript m) corresponding to the mid- interval pressure p. 
Thus , 
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Performing the integrations. 



RT 



AV 2 = - 



'm 



1 + r. 



m 



a W 
m g i 



log p + W.l 
e p - Ap/2 



m 



r Ap 
m 

b pi 
m m 



( 17 ) 



Equation (17) is the final form of the momentum equation. 

3 . Energy 

The enthalpy change of the mixture between state 1 

(the beginning of pressure interval Ap) and state 2 (the 

end of the interval) can be evaluated in two steps: 

(1) phase change at p. , T , T t and (2) change to P~ , T t 

1 "l "i z L 2 

Tg 2 , at fixed composition. 

The enthalpy change for step 1 is 



AH 



1 



[ amount of A i 
vaporized 



» « 

amount of B 
vaporized 



amount of A 
and B 
vaporized 

m 



enthalpy required 

to vaporize and 

heat unit mass of 

A from T t to T T 

L gL J 

enthalpy required" 
to vaporize and 
heat unit mass of 
B from V T to T 
L i 



X 



’ kinetic energy 
required to 
accelerate unit 
mass from V to 
V L 1 

g l 
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or 



4B = (m - m ) CL + o (T - T )] 

g 2 g x i g 1 9 l l 



+ (m - % ) [L b- + c b (t - t ) ] 

g l 1 



+ (m - m ) (vj* - V? )/2 
g 2 g l g l X 1 



( 13 ) 



where L and c are latent heat and specific heat, respectively. 
Introducing more compact notation. 



AH, = Ain (L + c 5,T) 

1 a a, a 1 

g 1 g x 

Am 5 , V 2 

+ im b < L b, + c b 5 1 T > + ^2 

2 1 g x 



(19) 



The enthalpy change for step 2 is evaluated from the 
temperature, pressure, and velocity changes, with properties 
evaluated at mean T and p for the interval. 



AH 0 = m [c (T - T ) + i (V 2 - V 2 ) ] 



2 ?2 g m ?2 



2 g. 



+ m [c x (T - T ) + 

± 2 m 1 2 ± 1 p 



P 2 ~ Pi 



m 



+ j (V 2 - V 2 ) ] 

^ ± 2 L 1 



( 20 ) 



37 



( 21 ) 





+ 



AV 2 

g 

2 



+ m. 





3y Assumption 3, no work is done by the free-stream flow and 
no heat is transferred to it. Hence, 



AH 1 + AH 2 = 0 



( 22 ) 



Substituting Eqs. (19) and (21) into Eq. (22) and 
2S the 
‘ AV 2 



solving for AT^ gives the energy equation for the mixture: 



AT 

g c . 



2 “ + r 2 



3 m 



<x 

\ m pl m 



Am o.V 2 Alil a 

+ -,3- 1 - . + 



2 * 



^2 



-r—1 (L + c 6.T) 
m a , a 1 

g-, i g-, 



Am, 



+ -r— 2. L, + c, 6.T 
m b. b 1 

9? 1 



(23) 



4 . Drag 

Although no force other than pressure acts on the 
free-stream flow as a whole, a drag force exists between 
the phases. Hence, a second momentum equation must be 
written using as the control volume the boundary between 
the phases. 

The two forces acting on each liquid drop are the 
buoyancy due to the pressure gradient and the drag due to 
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the relative gas velocity. The sum of these is equal to 
the mass times the acceleration of the drop. Thus, for a 
single drop. 



[dynamic pressure of 


y 


» < 

drag co- 


y 


7 * 

frontal area 


[relative gas flow 

• 


a 

i 


efficient 


A. 


of drop 



[ volume x pressure 
of drop gradient 

i * 



mass 
of drop 



X 



accelera- 
tion of drop 



or 




C ttD 2 ttD 2 dp 
4 ~ 6 dx 





(24) 



The absolute value sign in the first term makes the 
drag force positive when and negative when V < V^. 

Solving Eq. (24) for aV L , 



dv i = 



SPg I s | sV 2 C p dx 
4P 1 V 1 D 



_ dp _. 

Pl v l 



(25) 



Differentiating Eq. (12) , dV^ can also be expressed in terms 
of s, r, and V. Thus, 



dV 



1 



bdV + V 




sdr 

+ FT 2 



1 



ds 1 
+ r i 



•4 



(26) 



Solving for ds , 



ds 



b(l + r) dV 
V 



sdr 
1 + r 



(1 + r)dV x 
V 



(27) 
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Substituting dV from Eq. (25) , noting that dV = dv2/2V, 

Li 

using Eq. (12) , and writing for a finite increment, results in, 



As = 



b (1 + r ) AV 2 (1 + rJAp s Ar 
m m + _ m ^ + m 



2V 2 

m 



b Pt 

irr 1 m 



1 + r 
m 



m 



3p s | s C n (1 + r ) Ax 
g ‘ m 1 m D m 

m m 

4b p, D 
m I 

m 



(28) 



This is the drag equation employed when x is specified 
as a function of p. 

Solving Eq. (28) for Ax yields the required alterna- 
tive equation: 



Ax = 



4D 



3p s s C n V 2 
g m m D 
1 m 



b 2 p. 

m^l 



AV 2 



Ap + 



m 



+ 



b p. v2 . 

ml m / s Ar 
m / m 

1 + r m .1 + r 
m ^ 



m 




(29) 



5 . Heat Transfer 

Although no heat is transferred to the mixture as 
a whole, heat transfer exists between the phases. Hence, 
a second energy equation must be written using as the 
control volume the boundary between the phases. 

The work dW done on the liquid is that due to drag 
by the gas. (Only work done by shear or shaft forces is 
included in dW when writing the First Law for a control 
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volume) . Multiplying Eq. (24) by the number flow rate of 
drops N = 6m L /TTD^p L , the drag force on that quantity of 
liquid is 



F , = 



N 



8 P g 



V 



V C_.ttD 
1 g D 



2 _ 



m. 



d£ 

dx 



m i v i 



^i 

dx 



(30) 



The work done on the liquid is 



-dW 



F d dx 



- m/SE + 

Ml 



av ^ 

2 ) 



(31) 



The heat dQ transferred from the liquid is made up 
of two parts: (1) the convective cooling due to the 

temperature difference between the liquid and gas and (2) 
the evaporative cooling due to the latent heat supplied to 
the liquid vaporized. The convective cooling is 



-dQ c = hA d N(T 1 - T )dt 



(32) 



2 

where h is the heat-transfer coefficient. A, = ttD is the 

d 

surface area of a drop, and dt = dx/V L is the time required 
to traverse dx. Thus, 



-dZ 



c 



6hm 1 (T g - 



T^) dx 



Dp 



1 V 1 



(33) 



The evaporative cooling is 



-dQ v = L a cta a + 

g 



g 



(34) 
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The change in enthalpy of the liquid over the pressure 



increment dp is 



dH ajh.fc.dT. + 

1 \ 1 1 P 1 2 , 



(35) 



Substituting Eqs. (31), (33), (34), and (35) into the 

steady-flow energy equation dQ - dW = dH, the result is 



6hm^5Tdx 

Dp 1 V 1 



- L dm - L, dm, = m n c, dT. 
a a bp 111 

9 9 



(36) 



where 5T = T - T_ . 

g L 

Writing for a finite interval, the final form of the 
drop heat-transfer equation is 

1 



Am = — 
1 c, 



m 



6h6 TAx 
m 



Am 



Am, 



Dpi 



V, 



- L 



3- Jh-i 
m -L 



- L 



in m m m — 

Equations (1), (17), (23), (28) and (37) are the five 

equations that must be solved simultaneously to obtain the 
values of the five dependent variables D, T , T^, v , and 
V L as a function of the independent variable p. To carry 
out the solution all quantities in the equations must be 
expressed in terms of these six variables. 

The preceding equations form the basis for the 
mathematical model which is used to predict, based on 
inlet conditions, the exit velocity, and temperature of the 
mixture. These equations also form the basis for the model 
which provides the optimum nozzle shape given a set of inlet 



b m 1 
m i 



(37) 
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conditions. Some additional relationships are, however, 
required. These are: 

1. Phase properties — to establish the mass ratio, 
mass flow rate ratio of gas to liquid, and the 
thermal conductivity of the mixture. 

2. Liquid drop drag coef f icients . 

3. The liquid drop heat transfer coefficients. 

4. Boundary layer momentum thickness and displace- 
ment thickness. 

5. Skin friction coefficient. 

These five additional relationships are developed in detail 
in Ref. 1. 

The preceding "basic" model in conjunction with the 
additional equations form the core for four computer programs 
which provide real and isentropic flow output parameters, 
and an optimized nozzle contour given a set of inlet 
conditions . 

Reference 1 reported that the optimum nozzle profile 
has an elongated throat region, that the required shape is 
insensitive to drop diameter and nozzle length, and that the 
exit velocity is insensitive to departures from the optimum 
shape. For given fluids and pressures, the exit velocity 
was found to be insensitive to nozzle length and flow rate 
beyond certain minimum values. This insensitivity of 
nozzle exit velocity to nozzle shape and size results from 
the inverse dependence of drop size on slip velocity, which 
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apparently acts to adjust the slip to a value that holds the 
exit velocity in the range 75 to 90% of isentropic. Thus, 
it is difficult to design either a very good or a very poor 
two-phase nozzle. 
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V. EXPERIMENTAL APPARATUS 



The experimental apparatus was designed to provide a 
facility to study the performance of a Bi-phase nozzle. 

The basic requirements were to allow overall nozzle per- 
formance evaluation as well as optical diagnostic probing 
of the Bi-phase flow-field. 

The experimental system can be conveniently grouped into 
five subsystems. These are: the nozzle^ proper , nozzle 

assembly, air supply system, liquid injection system, and 
nozzle exhaust plane instrumentation. Each subsystem is 
described in the following sections. Figure 8 is an overall 
system drawing. 

A. NOZZLE 

The nozzle was designed based on recommendations given 

in Refs. 1 and 5. It is convergent-divergent, 12 in. long 

2 2 

with a throat area of 0.25 in. . Inlet area is 2.0 in. . 

The throat is located 4 in. from theinlet (1/3 of the over- 
all length) . The nozzle is constructed by "sandwiching" 
two h inch thick machined aluminum nozzle profile plates 
between h inch plexiglass plates (Figs. 9 and 10) . The 
aluminum nozzle plates are located at the end of a 30 inch 
long "test" section and are easily replaceable. The test 
section provides the necessary space for pressure taps, 
liquid injection tubes, and other instrumentation and 
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Fig. 8. Overall System Drawing 
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Liquid Injection System 
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Fig. 10. Nozzle Assembly Photograph. 
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configuration options. The ability to replace the nozzle 
profile plates allows for modification of the nozzle 
configuration . 

B. AIR SYSTEM (Fig. 11) 

Two 117 ft. ^ tanks supplied with compressed air from an 
Ingersol Rand three cylinder air compressor provide the air 
storage volume required to support nozzle operation. The 
nozzle is supplied with air at variable pressure via 3” ID 
piping. The air supply piping contains two pressure control 
valves, a solenoid actuated nitrogen operated 3 in. ball 
valve, and a flow measurement orifice plate. The orifice 
flow coefficient was determined from Ref. 6, Table 4. The 
air mass flow rate can be determined from the equation 

W, = 359 K d 2 Fa Y /h y 
h w 1 

where 

- air mass flow rate in LBm/HR 

K - orifice flow coefficient Table 4 [Ref. 6] 
d - orifice diameter 

F - thermal expansion factor Fig. 38 [Ref. 6] 

cl 

y - expansion factor. Fig. 40b [Ref. 6] 

h - differential pressure in H_0 at 68°F 
w r 2 

y - Specific weight of the flowing fluid at the inlet 

3 

side of the primary element (LB/FT ) 
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Fig. 11. Gas Flow Drawing. 
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and 



, _ m m o 

w 62.317 

where 

h - inches of differential of inches of manometer fluid 
m 

Y m - specific weight of manometer fluid LB/FT^ 

Y q - specific weight of fluid separating manometer from 
flowing fluid 

h w - differential head, in. of 1^0 at 68°F 

And 

Y = Y a 
Where 

Y, = specific weight of dry air at actual temperature. 

cl 

For an example calculation, see Appendix 2, example 5 from 
Ref. 6. 

C. LIQUID INJECTION SYSTEM (Fig. 12) 

The liquid injection system consists of a pressure 
vessel to contain the water which is pressurized from a 
nitrogen tank. The pressure vessel supplies pressurized 
water to the injection tube via a solenoid stop valve, a 
flow control valve and a flowmeter. The flowmeter was 
calibrated using an accurate scale and a stop watch. 

The calibration data and chart are recorded in Fig. 13. 

The liquid injector is a 0.25 in. brass tube inserted in 
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Fig. 12. Liquid Flow Drawing. 
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1!0 VAC 






Flow Rate L3 m /min 



Fig. 13. Liquid Flowmeter Calibration Curve. 
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the 3 in. ID air supply pipe just upstream of the flange 
connection to the test section. The injector tube is 
drilled with 5 holes facing the test section entrance. 

It has been proposed that more efficient liquid atomization 
could improve nozzle efficiency; therefore, the drilled 
holes were made as small as possible consistent with 
achieving a significant liquid mass flow rate. 

D. NOZZLE EXHAUST INSTRUMENTATION 

A cylindrical drag body, which is attached to a rectan- 
gular cantilever beam, is placed in the nozzle exhaust 
flow stream. The drag body unit is attached to a bracket 
which is attached to a rectangular transversing plate 
(Figs. 14 and 15) . The transversing plate is screw driven 
to move in a transverse direction across the nozzle exhaust. 
This movement causes the cylinder to move through the 
entire nozzle exhaust flow stream. By placing strain gages 
on the drag body unit, the cantilever beam deflection can 
be determined and related to the force on the cylindrical 
portion of the drag body. That force, in turn, can be 
related to the local velocity of the flow stream and these 
velocities plotted to give a velocity distribution at the 
nozzle exhaust. A similar plot can be constructed for the 
momentum distribution 

u 2 

Drag (LB f ) = C D A £y- 
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Fig. 15. 



Instrumentation Photograph. 
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where 



C D - Drag coefficient of the cylinder 

- (Smooth cylinder 1/8 in. in diameter) 

A - Projected area of the cylinder 1 / 2 " X 1/8" = 1/16 in.^ 
p - Density of the mixture 

U - Local velocity of the exhaust flow stream 
Solving for U: 

TT [2 X Drag force 
u ^/ — ^ 

Figures 16 and 17 are calibration data for the drag cylinder. 



57 



Hewlett Packard 92/0 toot. 
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Fig. 16. Drag Cylinder Calibration Curve. 
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Fig. 17. Drag Cylinder Calibration Curve. 
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VI . CONCLUSIONS 



Successful Bi-phase nozzle operation was achieved and 
sustained under controlled conditions. Liquid injection, 
nozzle, and instrumentation system all operate satisfactori- 
ly. The liquid breakup and acceleration in the nozzle 
passage was observed. 

With the nozzle pressures employed, i.e., 20, 30, 40, 
and 50 psig an instability was observed in the downstream 
portion of the nozzle passage. It was apparent that, for 
the nozzle geometry installed, there was a transition from 
supersonic flow in the throat region to subsonic conditions 
in the diverging downstream passage. It is apparent that 
at this point the passage converts from a supersonic nozzle 
to a subsonic diffuser. The transition point appeared to 
be unsteady, initiating alternatively from the sides of the 
nozzle. The position at which transition takes place could 
be controlled at will. As the liquid mass ratio and/or 
upstream pressure were increased, the passage tended to 
perform as a supersonic nozzle, i.e., the supersonic- 
subsonic transition was driven further downstream towards 
the exit. This is as one would expect. As the liquid 
mass ratio is increased, the sonic velocity decreases and 
the passage tends to a supersonic nozzle. Also, as the 
pressure ratio increases, the shock or transition to subsonic 
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conditions is driven downstream. An exit plane survey 
(Fig. 18) using the instrumented drag cylinder confirms 
the existance of a transition point somewhere in the diverg- 
ing portion of the flow passage. Thus, one is led to 
conclude that the role of change of flow passage area is 
too high for the test pressure ratios. From Fig. 18 it 
is also evident that the flow tended to oscillate from one 
side of the nozzle to the other. This was confirmed by 
visual observations as well as audible pulsations. 

It would be desirable to develop a degree of flexibility 
in changing or altering the nozzle profiles, particularly 
the diverging passage. Originally it was envisioned to 
replace the nozzle inserts; however, it may be more expedi- 
tious to provide an infinitely and continuously variable 
system. This may be achieved by a flexure system controlled 
by jack screws in the exit plane. The second desirable 
aspect would be the capability to diagnostically probe the 
flow field in the passage proper. A laser anemometer may 
provide the capability; however, no information exists on 
the use of LDA in Bi-phase flow. 
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Fig. 18. Nozzle Exit Force Distributions. 
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